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Abstract 

We exhibit two distinct renormalization scenarios in many-parameter 
families of piecewise isometries (PWI) of a rhombus. The rotational 
component, defined over the quadratic field K = Q(\/5), is fixed. The 
translations are specified by affine functions of the parameters, with 
coefficients in K. In each case the parameters range over a convex 
domain. 

In one scenario the PWI is self-similar if and only if one parameter 
belongs to IK, while the other is free. Such a continuous self-similarity 
is due to the possibility of merging adjacent atoms of an induced PWI, 
a common phenomenon in the Rauzy-Veech induction for interval ex¬ 
change transformations. 

In the second scenario, the phase space splits into several disjoint 
(non-convex) invariant components. We show that each component 
has continuous self-similarity, but due to the transversality of the cor¬ 
responding foliations, full self-similarity in phase space is achieved if 
and only if both parameters belong to K. 

All our computations are exact, using algebraic numbers. 

August 25, 2015 
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1 Introduction 


This work represents a first study of renormalization in many-parameter 
families of planar piecewise isometries (PWI). These are maps of polygonal 
domains partitioned into convex sub-domains —called atoms— in such a 
way that the restriction of the map to each atom is an isometry. The first- 
return map to any convex sub-domain D is a new PWI, called the induced 
PWI onD. If by repeating the induction we obtain a sub-system conjugate 
to the original one via a suitable group of isometries and homotheties, then 
we consider the original PWI to be renormalizable. 

Recent work on renormalization in two-dimensional parametric families 
concerned one-parameter deformations of the translational part of a PWI. 

110,16,23 . Induction is accompanied by a transformation s i —> r(s), where s 


is the parameter and r is the renormalization function. Self-similarity then 
corresponds to the periodic points of r. 

This setting is analogous to Rauzy-Veech induction for interval exchange 
transformations (IET’s, see |21,24, 26j), which are one-dimensional PWI’s. 
For an IET, the parameters are a vector of sub-interval lengths together 
with a permutation, and the renormalization acts on the lengths via an 
integral matrix. The fixed-point condition for self-similarity is an eigenvalue 
condition for a product of matrices. 

An arithmetical characterisation of self-similarity is provided by the 
Boshernitzan-Carroll theorem [5], which states that if an IET is defined 
over a quadratic number field (meaning that all intervals’ lengths belong 
to that field), then inducing on atoms results in only finitely many distinct 
IETs, up to scaling. In the case of two intervals (a rotation), this theorem 
reduces to Lagrange’s theorem on the eventual periodicity of the continued 
fractions coefficients of quadratic irrationals. However, unlike for continued 
fractions, there are self-similar IETs over fields of larger degree [3|. It is also 
known that in a uniquely ergodic self-similar IET, the scaling constant is a 


unit in a distinguished ring of algebraic integers 20 


In two dimensions general results are scarce [18, 19], All early results on 
renormalization concerned specific models of PWI’s defined over quadratic 
fields (the field of a PWI is determined by the entries of the rotation matrices 
and the translation vectors defining the isometries) 01 12||14[[22] . A more 
intricate form of renormalization has been found in a handful of cubic cases 
|9,[l5; 


The first results on parametric families concerned polygon-exchange trans¬ 
formations, due to Hooper [lO | (on the measure of the periodic and aperiodic 
sets in a two-parameter family of rectangle-exchange transformations) and 
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Schwartz |23| (on the renormalization group of a one-parameter family of 
polygon-exchange transformations). Subsequently, the present authors |l6 
studied two one-parameter families of piecewise isometries. Each family 
has a fixed rotational component defined over a quadratic field (Q(\/5) and 
Q(\/8), respectively), and parameter-dependent translations. It was shown 
that self-similarity occurs if and only if the parameter belongs to the respec¬ 
tive held. 

The mapping for the pentagonal model is shown in figure [lj For the 
parameter s restricted to a suitable interval, there is an induced PWI on a 
triangular sub-domain (the so-called base triangle ), which reproduces itself 
after scaling and the reparametrisation s i —> r(s). After an affine change of 
parameter, the function r was found to be of Liiroth type —a piecewise affine 
version of Gauss’s map |4,8 JTt] . In the pentagonal model, the discontinuities 
of r accumulate at the origin; in the octagonal case one has r = fof , where / 
has two accumulation points of discontinuities. In both cases r is expanding 
and preserves the Lebesgue measure. 



Figure 1: One-parameter rhombus map of the pentagonal model of 16 . The 


rhombus and all directions are fixed. The parameter shifts the inner boundaries of 
the atoms in such a way that all three atoms retain a reflection symmetry. 


The present work deals with two-parameter extensions of the one-parameter 
pentagonal model of [l6| (figure [2]). The roles of the parameters is made 
explicit in figure [3j where we employ a new co-ordinate system in order 
to restrict the arithmetic to the quadratic field K = Q(\/5), rather than 
a bi-quadratic field —see section 2.6 In the new co-ordinates, the atom 
boundary lines are of the form 


u 


• X = bi 


h = b ifi + bi,isi + b i}2 S2, bij€<Q(V 5) 


where si, S 2 are the parameters. We use these co-ordinates for calculations, 
reserving the original co-ordinates for graphics. 


3 








Figure 2: Two-parameter rhombus map. 


In section [2] below, we develop the necessary geometrical constructs: 
tiles, dressed domains, isometries, etc. In section [3] we review the renormal¬ 
ization theory of the base triangle developed in 16], with two additions. We 
extend the isometries to the atoms’ boundaries (which will be needed to glue 
together atoms in section [4]), and we introduce an improved renormalization 
scheme whose renormalization function has finitely many singularities. 

The first renormalization scenario is established in section |4j There we 
make three successive inductions on triangular sub-domains, to obtain a 
PWI whose eight distinct return orbits tile the rhombus, apart from the 
orbits of finitely many periodic domains. We identify a convex polygon II 
in the t. s parameter space within which the induced map has only three 
distinct isometries. By merging neighbouring atoms which share the same 
isometry, we recover the one-parameter PWI discussed in the previous sec¬ 
tion. As a result, self-similar dynamics occurs if and only if s E K, and since 
t is arbitrary, the self-similarity constraint corresponds to a foliation of II. 

The appearance of a hidden reduction of the number of atoms, which 
reveals itself only after induction, is present in a simpler form in the Rauzy- 
Veech induction of interval exchange transformations. In section [5] we show 
that the occurrence of free parameters in renormalizable IETs can be under¬ 
stood in terms of the properties of the associated translation surfaces [ll . 
In particular, there are self-similar IETs without free parameters, as long as 
the number of intervals is greater than three. 
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Figure 3: The two-parameter map conjugate to the rhombus map of figure [2| 

In view of this analogy, we conducted an extensive search for a non¬ 
degenerate two-paranreter family, using two-dimensional sections of a three- 
parameter system of the 27r/5 rhombus (sectionjb]). We found a weak form of 
non-degeneracy, resulting from the co-existence of two continuously renor- 
malizable families with transversally intersecting foliations. Each family 
corresponds to a (non-convex) invariant set in phase space, and hence to an 
ergodic component of the exceptional set. 

In many-parameter families of piecewise isometries over quadratic fields, 
the existence of an irreducible component in phase space which is self-similar 
without free parameters has yet to be established. 

Acknowledgements: JHL and FV would like to thank, respectively, the 
School of Mathematical Sciences at Queen Mary, University of London, and 
the Department of Physics of New York University, for their hospitality. 

2 Preliminaries 

Throughout this paper, we let 

a = y/b, co = (a + l)/2, (3 = w _1 = (cc — l)/2. (1) 

The arithmetical environment is the quadratic field Q(u;) with its ring of 
integers Z[u], given by 

Q(w) = {x + ya : x, y e Q}, Z[u\ = {m + nui : m,n € Z}. (2) 
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The number u>, which is the fundamental unit in Z[w] (see J6, chapter 6]), 
will determine the scaling under renormalization. The number /3 = u — 1 is 
also a unit. 

2.1 Planar objects 


A tile X with n edges is 

a convex polygon defined by the half-plane condi- 

tions 


u m , • x < 

bi (excluded edge) 

or 

z = 1,..., n, (3) 

u m , • X ^ 

bi (included edge) 


where x = (x,y), bi G M, and the u m are the vectors 


( 2irm 2-Km\ 

u m = I cos —— ,sm —J me{ 0, ...,4j. (4) 

For the ith edge, defined by u mi ■ x = bi, we introduce an index e,, where 
e* = —1 if the edge is included in X , and e* = 1 if it is excluded. We then 
represent X as a triple of n- vectors 

X = [(mi,... ,m n ), (ei,... ,e n ), (h,.. .,b n )}. (5) 

We shall assume that n is minimal, namely that X is not definable by fewer 
conditions. 

A tiling X is a set of disjoint tiles, 

X = {Xi ,..., X N } 

and is associated with a domain X (union of tiles) 

N 

x=\Jx k . 

k= 1 

Note that a domain need not be convex, or even connected. Note further that 
thanks to if a pair of tiles have disjoint interiors but share a common 
boundary segment, that segment belongs to one and only one tile of the 
pair. This allows the possibility of gluing together adjacent tiles without 
disturbing the inclusion relation of the respective edges. 
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2.2 Similarity group 

The transformation properties of planar objects are provided by a group 05 
which comprises the rotations and reflections of the symmetry group of the 
regular pentagon (the dihedral group D§) together with translations in K 2 
and real scale transformations. 

We adopt the following notation: 

U m : reflection about the line generated by u m . 

R m : rotation by the angle 2rmr/5. 

Tcj: translation by d e K 2 . 

S^: scaling by 77 E M+. 

We write X ~ y to indicate that X is similar to y, i.e., that X = G(Y) for 
some G 6 0. As 25 is a group, this is an equivalence relation. Within & we 
distinguish two important subgroups: the isometry group 3 generated by ro¬ 
tations, reflections, and translations, and the dynamical group 3+, generated 
by rotations and translations. 

2.3 Dressed domains and sub-domains 

A dressed domain is a triple 


Y = (X,X,p), (6) 

where X = {X \,..., Xn} is a tiling of the domain X, and p = {pi,..., pat}, 
where p* E 3+ is an orientation-preserving isometry acting on the tile X^. 
Under the action of G E <5, a dressed domain X transforms as 

G(A) = G(X, X, p) = (G(X), {GprO,..., G(X k )},G o p o G” 1 ) 

where the conjugacy acts componentwise. To emphasize the association of 
a mapping p with a particular dressed domain X, we use the notation p x . 

Let X = (X, X, p x ) be a dressed domain, and let Y be a sub-domain of 
X. We denote by Py the first-return map on Y induced by p x . We call the 
resulting dressed domain y = (Y, Y , Py) a dressed sub-domain of X , and 
write 

y <x. ( 7 ) 

The dressed sub-domain relation 0 is scale invariant , namely invariant 
under an homothety. Indeed, if denotes scaling by a factor 77 , then in the 
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data ([5]) specifying a tile, the orientations mk remain unchanged, while the 
pentagonal coordinates bk scale by 77 . Moreover, the identity 

S^TfjRn = T^dRjjS^. 

shows that the piecewise isometries p scale in the same way. We conclude 
that the subdomain relation Q is preserved if the dressed domain parame¬ 
ters are scaled by the same factor for both members. 


2.4 Parametric dressed domains 

In this article we consider continuously deformable dressed domains X = 
X(s) called parametric dressed domains , depending on a real parameter 
vector s = (s\,..., s p ). 

These are domains whose tiles and image tiles Pk{Xk) depend on s 
only via the coefficients bi, while the parameters n, mi and e* remain fixed 
[see ([5])]. We shall require that the bi s be affine functions of si,..., s p , with 
coefficients in Q(w). Algebraically, this is expressed as 


bi G S S — Q(w) + Q(w)si + • • • + Q(w)sp, (8) 

where si,..., s p are regarded as indeterminates. The set § is is a (p + 1)- 
dimensional vector space over Q(uj) (a Q(o;)-module). 

The condition ([8]) gives us affine functions bi : —> M 

bi(s\, ..., s p ) = bifi + fejqsi + • • • + bi y pSp bij G Q(w). (9) 


We define the bifurcation-free set n(T) to be the maximal open set such 
that all of the edges of all Xk(s) have non-zero lengths. Note that other 
types of bifurcations may occur if X is embedded within a larger domain 
(see section 2 . 6 .) 


2.5 Renormalizable dressed domains 

A parametric dressed domain X (s) is said to be renormalizable over an open 
domain II C M p if there exists a piecewise smooth map r : II —> II such that 
for every choice of s G r^ 1 (II) the dressed domain X(s) has a dressed subdo¬ 
main y similar to A(r(s)) which satisfies the recursive tiling property. The 
function r depends only on s, a requirement of scale invariance. In general, 
we have y = Ti( s )('S), where i is a discrete index. The set r _ 1 (II) need 
not be connected (even if i is constant), each connected component being a 
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bifurcation-free domain of T- (To extend the renormalization function r to 
the closure of n, one must include bifurcation parameter values, as in |jl6|.) 

If s = so is eventually periodic under r, then we say that T’(so) is self¬ 
similar. A self-similar system has an induced sub-system which reproduces 
itself on a smaller scale under induction. 

Let a parametric dressed domain X(s) have induced Xj(s), such that, 
for j = 1,... ,N we have: (z) Xj is renormalizable over a domain IP,; (zz) 
the Xj recursively tile X; (in) the II j have non-empty intersection II. Then 
we still consider X renormalizable over II. 

The definition of renormalizability given above is tailored to our model; it 
is not the most general possible, and it is local in parameter space. We allow 
T to depend on a discrete index (as in Rauzy induction for interval-exchange 
transformations —see section[5]) to obtain a simpler renormalization function 
r (section [3]). We only require X to be eventually renormalizable, and we 
allow X to have sub-domains with independent renormalization schemes 
(which is a common phenomenon, see section [6]). 


2.6 Computations 

For computations, we use the Mathematica®procedures listed in the Elec¬ 
tronic Supplement [7j. All computations reported in this work are exact, 
employing integer and polynomial arithmetic, and the symbolic representa¬ 
tion of algebraic numbers. 


The geometrical objects defined in section 2.1 require arithmetic in a bi¬ 
quadratic field, since only the first component of the vectors u m is in Q(w). 
To circumvent this difficulty, we conjugate our PWI to a map of a square 
where the clockwise rotation 2tt/5 is represented by the following matrix 
over Z[w\ 

0 1 
-1 p 

where (5 was defined in 0- (This is still a PWI with respect to a non- 
Euclidean metric.) In the new co-ordinates, the vectors u m become 


{( 1 , 0 ), ( 0 , 1 ), (- 1 , 0 ), (- 0 ,- 0 ), ( 0 ,- 1 )} 

which belong to Z[cu] 2 . With this representation, all of our calculations can 
be performed within the module S defined in ([8]). We shall still display 
our figures in the original coordinates, where geometric relations (especially 
reflection symmetries) are more apparent to the eye. 

In constructing a return map orbit of a domain X(s) by direct iteration, 
one determines inclusion and disjointness relations among domains, which 
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require evaluations of inequalities <©■ Since the latter are expressed by 
affine functions of the parameter s in some polytope n, it suffices to check 
the inequalities on the boundary of II. All these boundary values belong to 
the field Q(w), and the inequalities can be reduced to integer inequalities. 

Typically, X will be immersed in a larger domain y (an atom, say). 
Therefore, in addition to the intrinsic bifurcation-free polytope n(T) defined 
in section 2.4, one must also consider the polytope II (X,y) defined by the 
inclusion X(s) C T(s), as well as intersection of these polytopes. 

The recursive tiling property defined in section 2.5 is established by 
adding up the areas of the tiles of all the orbits, and comparing it with the 
total area of the parent domain. 

With these techniques, we are able to establish rigorously a variety of 
statements valid over convex sets in parameter space. 


3 Base triangle 


The base triangle is the simplest one-parameter renormalizable piecewise 
isometry associated with rotations by 27 t/ 5; it is self-similar precisely for 
parameter in the quadratic field Q(\/5). It was instrumental to the proof 
of renormalizability of a one-parameter rhombus map in 
appear again in the many-parameter versions presented here 
We develop a variant of the model presented in 


16 , and it will 


16 , which includes 


boundary segments in the tiles and an improved renormalization scheme. 
The base triangle B prototype is the following dressed domain (see figure 
0 : 
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B — ( B, (Bi, B 2l B 3 ), ( pi,P2,P3 )) 


where 

B = [(1,0, 2), (-T, 1,1), (t — w 2 ,0,0)], 

Bi = [(0, 2,3), (1,1,1), (0,0, w -ut)\, (10) 

B 2 = [(1,4, 0, 3, 2), (—1, 1,1, —1, 1), (r — a; 2 , cj 2 — lot , 0, w — lot , 0)], 

^3 = [(1,0,4), (—1,1, —1), (r — w 2 ,0,u; 2 — ur)]. 

The dynamics is given by a local reflection of each atom about its own 
symmetry axis, followed by a global reflection about the symmetry axis of 
B, which can be written as: 


Pi 



P 2 

= T (0,o; 2 t— a; 3 ) R 3 

(11) 

P 3 

f(a;r— 2uj,lu 2 t— 2u> 2 ) ^2- 



Here we have chosen a coordinate system such that the peak of the isosceles 
triangle is at the origin and the altitude of the atom B3 is the parameter 
r, which varies over the interval (0,1) without the occurrence of a bifur¬ 
cation. This parameter (together with time-reversal invariance) determines 
the scale-invariant properties of the dressed domain, since it is related to 
the ratio 77 of altitudes of B 3 and B by the formula 

T 

V = ~2 -• 

(jJ z — T 

As t varies from 0 to 1, r] increases from 0 to f3. 

The edges of the domain B are included or excluded as stipulated in 
section [2| a vertex joining two included edges is included, but is excluded 
otherwise. The renormalizability analysis will also require a second base 
triangle B , differing from B by a change of sign of all edge coordinates and 
translation vectors, as well as of the respective ej. The dressed domains B 
and B are ©-inequivalent: not only do they have different boundary con¬ 
ditions, but their interiors differ by a rotation by 7 r, not an element of the 
similarity group. 

The renormalizability analysis for the base triangle is summarized in the 
following lemma: 

Lemma 1 Let B be as above. The following holds: 

(i) For 0 < r < f3 2 , B has a dressed subdomain B\ ~ B which is scaled by 
a factor (1 — t)/(uj 2 — r) and has shape parameter r(r) = w 2 r. 
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(ii) For f3 2 < t < (3, B has a dressed subdomain £>2 ~ £> which is scaled by 
a factor r/(w 2 — r) and has shape parameter r(r) = u 3 (/3 — r). 

(Hi) For (5 < t < 1, B has a dressed subdomain £>3 ~ B which is scaled by 
a factor t/{uj 2 — r) and has shape parameter r(r) = w 2 (l — r). 



Figure 5: Renormalization function r(r) for base triangles. 


The renormalization function r has three branches (see hgure[5]). In cases 
(ii) and (iii) one induces on the atom S 3 , over two disjoint bifurcation-free 
parameter ranges. Since the size of S 3 vanishes as r approaches 0, in the 
range (i) we induce on the triangle [( 1 , 0 , 2 ), (— 1 , 1 , 1 ), (/3 2 r — 1 , 0 , 0 )], which 
is not an atom. This device prevents the occurrence of infinitely many 
singularities in the renormalization function found in [16;. 

In each case, the return orbits of the atoms of £>,;, together with a finite 
number of periodic tiles, completely tile the triangle S. For B , the pre¬ 
scriptions (i)-(iii) hold with the roles of B and B exchanged. The induction 
relations are represented as the graph in figure [6j 

As in 16 , the proof of Lemma 1 is by direct iteration, as discussed in 


section 2.6 The main difference in the two computational algorithms lies 


in the procedures used to verify inclusion and disjointness of tiles (see 0 )- 
Specifically, checking the sub-polygon relation X C Y requires verifying that 
no included vertex of X has landed on an excluded edge of Y. Similarly, to 
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Figure 6: Renormalization graph for base triangles. A directed link from X to 
Y indicates that X has an induced dressed subdomain equivalent to Y, subject to 
the the indicated parameter constraint. 


decide that X and Y are disjoint, one must check that no included vertex 
of either tile lies on an included edge of the other. 

With reference to lemma [lj we remark that the base triange B, with 
its definition extended to the two-atom limiting cases r = 0,1, is in fact 
renormalizable also at the parameter values 0,/3 2 ,/3, 1, with r(r) = 0 in all 
these cases (see [7] for the calculations). These additional parameter values 
are needed to make B renormalizable over the whole interval [0,1] n Q(w). 
We shall use this property in sections [4] and [6] 


4 Continuous self-similarity 


We now turn to the two-parameter rhombus map introduced in section [l] 
—see figures [2] and [3j In suitable coordinates, the dressed domain is given 
by 

1Z = (R, (Ri, ... ,i? 5 ), (pr 1 ,. .., Pr 5 )), 

with (see figure [2]) 


R = [(0,1,0,1), (—1, —1,1,1), (—t, —s, 1 — i, 1 — s)], 

Ri = [(0,2,1), (—1, —1,1), (—f, —s, 1 — s)], 

R 2 = [(0,1, 2,0,1, 2), (—1, —1, —1,1,1,1), (—t, —t, —1 — s, 1 — t, 1 — s, —s)], 

R 3 = [(0, 2,1), (1,1, —1), (1 — t,—l — s, —i)], (12) 

R 4 = [(0,1, 2,1), (—1, —1, —1,1), (—t, —s, —1 — s, —t)], 

R 5 = [(1,0,1, 2), (—1,1,1,1), (—s, 1 — t, —t, — 1 — s)] 


PRi = T (0 ,o) R 4, PR 2 = T(o,l) R 4, PR 3 = 1(0,2) R 4, (13) 

PRa = 1(1,1) R l- PR 5 = 1(1,2) R 4- 

The corresponding bifurcation-free parametric domain 11(77), defined in 
section 2.4 is found to be the triangle with vertices at (0, 0), (—1/a, —1/a), 
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and (/3/a, —/3/a). On the boundary of 11(7?.) given by with s = t, the dressed 
domain 7? collapses into the one-parameter pentagonal model of [l6j, and 
hence is self-similar for all s E Q(w) within a suitable interval. 



Figure 7: The first two steps of the triple induction 7? > 7?3 > 7?3i > A. The third 
step produces the dressed domain A shown in figure [8] 

Our goal is to determine a parametric dressed domain A with bifurcation- 
free subdomain 11(^4) C 11(7?) over which 7? is renormalizable. To this end, 
we choose a specific parameter pair close to the s = t boundary: (so, to) = 
(—19/200, —1/10), and we establish that at this value the renomalization is 
amenable to a complete analysis (with computer assistance). Specifically, 
we consider a three-step induction, first on the triangular atom R 8 , followed 
by two inductions on sub-triangles, as shown in figure [7} The last induction 
produces the dressed domain A, shown in figure [8j which is given by: 

A = (A, (Ai,.. .,A 8 ), (pa 1} - ■ ■ 1 PAg))i (14) 

with 

A = [(2,1,0), (1, —1,1), (—1 — s, /3 4 — s, 1 — s)] 

Ai = [(1,0,4), (—1,1, —1), (/3 4 — t, 1 — s, a/3 4 — t)], 

A 2 = [(1,4,0,4), (—1,1,1, —1), (/3 4 — t, a/3 4 — t,l — s, a/3 4 — s)], 

A 3 = [(0,1,4,1), (1,1, —1, —1), (1 — s, /3 4 — t, a/3 4 — t,/3 4 — s)], 

A 4 = [(4,1,4,1), (1,1,-1,-1), (a/3 4 -t,a/3 4 - s,/3 4 - s)\, 

A 5 = [(2,1, 4,0, 3), (1,-1,1,1,-1) (15) 

(—1 — s, f3 A — t, a/3 4 — s, 1 — s, —4/3 3 — t)], 
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= [(2,3,0,3), (1,1,1, —1), (—1 — s, — 4/3 3 — t, l — s, —4/3 3 — s)], 

A 7 = [(2,1,4,1), (1, —1,1,1), (—1 — s, /3 4 — s, a/3 4 — s,/3 4 — t)], 

A 8 = [(2,3,0), (1,1,1), (—1 — s, —4/3 3 — s, 1 — s)], 

PA\ = PA 2 = PA 3 = PA± = T(8£3 ,_2/3 s ) R 2, 

PA 3 PA 6 PA-j T(2,l+/3 5 ) R 3) (16) 

Pa 8 = T( 2 _ ( a6 i _ /3 5) R 2 . 




Figure 8: The dressed domain A, with its 8 atoms numbered as in (15). The 
boundaries of the composite atoms C\,C 2 ,C^ are coloured red, green, and blue, 
respectively. 


We find that 11(^4) is the triangle with vertices 


(-l+ 2 ,-l+ 2 ), (*(ll-5a), *(-25 + lla)), (*(11 - 5a), *(11 - 5a)), 

a a 2 4 2 2 

shown in figure [9j One verifies that 11(A) is adjacent to the line s = t and 
that ( sq, to ) lies in its interior. 

Using direct iteration, we verify that for all (s, t ) 6 11(A) the return 
orbits of the eight atoms of A, together with those of 13 periodic tiles, 
completely tile the rhombus R (see figure 10). 


A decisive simplification of the analysis results from the observation that 
the four atoms of A are mapped by the same isometry, and 

hence, with regard to the first-return map to A, can be merged into a single 
triangular tile, A 1234 . Similarly, atoms A 5 , Aq, Aj can be merged into a single 
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-0.100-0.095-0.090-0.085 


Figure 9: The parametric domain n(A). 

reflection-symmetric pentagon, A 567 . The mergers have been suggested in 
the shading of the tiles in figure [ 8 } The dressed domain thus simplifies into 

C = (C,(C\, C r 2 , C y 3 ), (pc'i, PC 2 ,PC 3 )) 

def 

= (A (^-8,^-567,^1234), (PA S , PA a , PA t )) (17) 

Moreover, one verifies that over n(C) = n(A), we have C ~ B. The intrinsic 
shape parameter of C can be calculated from the ratio rjc of the altitude of 
C 3 to that of C: 

T c = = u 7 (as + ff 3 ). (18) 

1 + VC 

As we transverse 11(C) from left to right, s increases from 2/a — 1 to (11 — 
5a)/2, with tq increasing from 0 to 1. 

The issue of recursive tiling is now rather subtle. The rhombus is cer¬ 
tainly tiled by the return orbits of C\, C 2 , C3, and the 13 periodic tiles which 
arose in the induction on A. (see figure |T7)j) . However, the return paths are 
not the same for all tiles. In (say) C3, the tiles Ai, A 2 , A 3 , A 4 have four 
distinct 78-step return paths, which go their separate ways, but recombine 
eventually to form an atom of the dressed domain with a unique isometry. 
The coincidence of the return times is not necessary to the recombination, 
as these times could differ by any integer multiples of 5. As a result, the 
partition of C 3 into A\, A 2 , A 3 , A 4 is relevant to the recursive tiling of the 
original rhombus, but not to dynamical self-similarity. (We shall encounter 
again the same phenomenon —recombination with different return times— 
in section §) 
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Figure 10: Tiling of 1Z 3 by return orbits of the 8 atoms of A (coloured) and 7 
periodic tiles (grey). Note that Ai, A 2 , A 3 , A 4 , which comprise Ci, have distinct 
return orbits. 


The parameter pairs (s, t) G 11(C) corresponding to self-similarity for the 
rhombus map 1 Z are now determined by the self-similarity of the induced 
dressed subdomain C. In turn, the latter are the values of s for which the base 
triangle is self-similar, namely t^(s) G (0, l)nQ(o;), while t is unconstrained. 
Thus, by (18), 1Z is renormalizable in n(C) if and only if 

(s, t) G n(C) n (Q(w) x M). 


This is the main result of this section. 
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5 Continuous self-similarity in Rauzy induction 


Recombination of atoms, and the resulting appearance of a free parameter 
in self-similarity may seem a coincidental feature of planar PWI’s. This phe¬ 
nomenon is in fact common in the Rauzy-Veech analysis of renormalizable 
interval exchange transformations (IET’s) 121,24,25 . 


We fix a half-open interval O = [0, l ) and a partition of O into n half-open 
sub intervals flj. An IET is a piecewise isometry of 0 which is a transla¬ 
tion on each Oj. We represent it as a pair (n, A), where A = (Ai,..., A n ) 
is the vector of the lengths of the sub-intervals and n is the permuta¬ 
tion of {1,..., n} such that the intervals in the image appear in the order 
tt( 1), ... ,vr(n). 

We assume that it is irreducible in the sense that {1,..., k} is mapped 
into itself only if k = n. If we fix n, then ( 7 r, A) is a parametric PWI, 
with discrete and continuous parameters ir and A, respectively. IET’s which 
differ only by an overall translation or scale transformation are considered 
equivalent. 

The Rauzy-Veech induction on ( 77 , A) consists of inducing on the larger of 
the two intervals f^ 0) = [0, l — A n ) and = [0, l — X^-i, denoted by type 
0 and type 1 induction, respectively (the case is excluded from 

consideration, as in this case the map is not minimal). Induction corresponds 
to a map (7r, A) i->- ( 7 P, A'). Letting 


(cLi(ir), Ai(n) 1 A) = (7t',A') * = 0,1 

one finds that Aj(7r) _1 is an n x n integral matrix (see 120 for explicit 
expressions for m and Ai(ir)). 

The permutations 7r of n symbols are then represented as the vertices 
of the Rauzy graph. Each vertex has two outgoing and two incoming edges, 
associated with ai and aj 1 , respectively, for i = 0,1. The Rauzy classes 



Figure 11: Rauzy graph for n = 3. All permutations are degenerate, and their 
translation surface is a torus. Any renormalizable IET with three intervals will 
have a free parameter. 
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are the connected components of the graph. These IET’s are linked by 
a sequence of Rauzy inductions, and a self-similar IET corresponds to a 
path on a Rauzy graph which terminates in a closed circuit, 7Ti, 7r2 ,..., ir p . 
Transversing such a circuit produces an induced IET which is a rescaled 
version of the original one. Its length vector is an eigenvector of a product 
of matrices A lk (n}^ 1 . k = 1 ,... ,p, with a scale factor given by the corre¬ 
sponding eigenvalue. In figures 11 and 12 we display the Rauzy graphs for 
n = 3 and n = 4 125l section 6]. 






Figure 12: Rauzy graph for n = 4. All permutations of the upper component are 
degenerate. No permutation of the lower component has this property. 

Two of the three irreducible permutations of the Rauzy graph for n = 3 
(figure 0 and all of the permutations for the first class for n = 4 (hg- 
ure |12[) contain a consecutive pair: (... ,j,j + 1,...). For those IET’s, the 
consecutive intervals f lj and flj+i have the same translation vector, and 
hence the interval Qj U flj+i may be merged into a single interval of length 
A '■ = A j + A J+ i. The n-interval IET is thus equivalent to an n — 1-interval 
IET. If the latter is self-similar, then the original PWI is also self-similar for 
any choice of the parameter A j £ [0, A']. 

Accordingly, we say that a permutation is degenerate if it has consecutive 
pairs or if it acquires this property after a single induction. In the latter 
case, consecutive atoms of the child IET have distinct return paths in the 
tiling of the parent. 
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Such a degeneracy is best understood by representing an IET as Poincare 
section of a flow on a translation surface nm 24j]25|. The latter is a poly¬ 


gon with 2 n sides (n is the number of intervals), labelled according to the 
ordering of the intervals before and after the permutation. The sides that 
correspond to the same interval have equal length and are parallel, and they 
are to be identified. A rectilinear flow on the plane will develop conical 
singularities on the surface, in correspondence to the vertices of the 2 n- 
gon. While the translation surface is not unique, its genus and singularities 
(given by the total angle 2 n(m + 1) at the identified vertices) depend only 
on the Rauzy class. The removable singularities (m = 0) correspond to de¬ 
generate permutations, and they signal the appearance of free parameters 
in renormalizability. Since the translation surface does not change under 
induction, these structures depend only on the Rauzy class to which the 
permutation belongs. Furthermore, for any n ^ 4 there are both degenerate 
and non-degenerate permutations, e.g., (n,..., n — 1,1) and (n,n — 1,..., 1), 
respectively. 

The permutations which acquire consecutive pairs after induction are 
uniquely of the form (..., n, k + 1 ,,k) for some k < n — 1, and all four 
of its neighbours in a Rauzy graph have consecutive pairs, thanks to the 
relations 


a 0 ((.-n, k + 1,..., k)) = (...,, k + 1, k + 2,..., k), 


-i 


((...n,k + 1,...,&)) = (-,n-l,ra, 

ai ((...n, k + 1 ,..., k)) = (..., n, k, k + 1 , ...), 
a^ 1 ((...n, k + 1,..., k)) = (..., n, ..., k, k + 1). 


The n = 3 class provides the simplest illustration of this phenomenon. 
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6 Weakly-discrete self-similarity 


The analogy with Rauzy induction suggests that there might exist two- 
parameter planar PWI’s which do not admit self-similarity with free param¬ 
eters, meaning that both parameters would be algebraically constrained. We 
shall exhibit a weak form of this property, resulting from the coexistence of 
two systems with continuous self-similarity. 
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Figure 13: Four-atom rhombus PWI. 



Our starting point is the four-atom, three-parameter PWI of the 27r/5 

The interpretation of the parameters s, t, u is 


rhombus shown in figure 13 


made clear in the conjugate system of figure p~4| where the rhombus appears 
as a unit square. (A similar strategy can be pursued for the five-atom family 
of figures [2] and [3j but we found that the four-atom family is somewhat easier 
to work with.) 

In convenient coordinates, the dressed domain is 


K = (R, (Ri,...,R 4 ),{pi,...,P4)), 


with 

R = [(0,1,0,1), (—1, —1,1,1), (— f, —s, 1 — t, 1 — s)], 
Ri = [(0,1,2,1), (—1, —1, —1,1), (— t, —t, —s, 1 — s)], 
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Figure 14: Dressed domain 1Z. Dependence of atoms on parameters s,t, u is 
shown. 


R 2 = [(0,1, 2,1), (1,1,1, —1), (1 — i, 1 — s, —s, —t)], 

R 3 = [(0,1,2,1), (-1,-1, -l,l),(-t,-s,-l-s + u,-t)}, (19) 

R 4 = [(0, 1,2, 1), (1,1,1, —1), (1 — t, — t, —1 — s + u, — s)] 


Pi — T(0,0) R 4, 
P3 = 1(1,1—u) R 4 


P'2 = T(0,1) R 4, 

PA = 1(1,2 -u) R 4 • 


( 20 ) 


From figure 14, we see that the bifurcation-free domain n(7£) C is the 
polytope bounded by the planes s — t = 0, s — t = 1, u = 1, s + (3 2 t = 0, 
u — s — /3 2 t = 0, /3 2 — ( 3 2 s — t = 0, and 1 + f3 2 s + t — u = 0. 

This system has a simple one-parameter subsystem on the line L defined 
by s — t = fj 2 , u = (3. We shall consider a two-parameter perturbation of this 
subsystem in the plane u = j3 which intersects L. (We have also considered 
other planes, obtaining other manageable examples: see remarks at the end 
of this section.) 

Setting u = j3, the parameter polytope reduces to the hexagonal domain 
shown in figure[T t| (left). As done in section |4j we choose a parameter pair 
close to L lying within such a domain: (sq, to) = (2/5,1/25). By inducing on 
the trapezoidal atom R \, we obtain the parametric dressed domain T shown 
in figure [15} One readily verifies that the return orbits of the eight atoms 
of T completely tile R, so that the renormalizability of 1Z will follow from 
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Figure 15: Induced dressed domain J. 


that of T. We find that the complete tiling of F by renormalizable dressed 
sub-domains, given in figure [16J requires the return orbits of three dressed 
triangles F 2 , F 3 , plus seven periodic tiles Vi (five regular pentagons, one 
trapezoid, and one rhombus). 

Letting 

3 7 

n*(.F) = p|n( 7 i) (~)n(v t ) (21) 

1=1 1=1 

we find that II* (F) is the quadrilateral with vertices 

(/? 2 ,0), (/3 2 + /? 6 /«,/? 6 /«)), (/3 2 + £ 6 /a,/5 4 /«), (/3 2 , £ 4 /«) , 

shown in figure [Tt] (right). Note that one of the bounding edges of the 
parameter domain coincides with the line L which was the starting point of 
our perturbative exploration. 

The dressed domains T\ and Ti are base triangles equivalent to the 
prototype B , with respective shape parameters 

7~i = w 6 a(s — f3 2 ), 72 = uj A at. 

Examination of J-$ shows that its atoms with four and five sides share the 
same isometry (in spite of having different return paths, and even different 
return times on the rhombus), and hence can be merged for the purpose of 
testing renormalizability. After the merger, J -3 is also equivalent to B , with 
shape parameter 

73 = w 6 a(s - /3 2 ). 
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Figure 16: Tiling of F by return orbits of F\ (red), (green), F 3 (blue), and 
seven periodic tiles (grey). 


Since each T) is renormalizable for r* e Q(\/5), we conclude that T 
(hence 1Z) is renormalizable when all three shape parameters are in Q(\/5), 
i.e., when (s,t) is constrained to belong to II*(J 7 ) flQ(\/5) 2 . 

Each of the three renormalizable dressed domains Jy provides a sequence 
of nested coverings of a distinct invariant component of the exceptional 
set complementary to all periodic points of the rhombus. The number of 
distinct ergodic components of the exceptional set is thus at least three. For 
the model of section |4j on the other hand, we believe that there is a single 
ergodic component. 

In closing, we summarize briefly the results of our explorations of the 
three-parameter space of the four- or five-atom rhombus maps. 

Manageable renormalizations are likely to be found in systems specified 
by parameters of small height. [The height H{Q of the algebraic number ( = 
{m/n)+{m!/n')oj is defined as H( £) = max(|m|, |n|, |m'|, |n'|)]. Such was the 
case for the domain 1Z, and the two-parameter restriction u = 0 described 
above. We have considered other planes of small height: s — t — uju + /? = 0, 
s — t + u— 1 = 0, etc. Within such planes, and for carefully chosen parameter 
patches, we have encountered a number of two-parameter renormalizable 
models. These are characterized by a decomposition of the rhombus into 
N disjoint dressed domains, each tiled by the return orbits of a single base 
triangle (provided we ignore the “decorations” produced by the common 


24 













Figure 17: Left: The bifurcation-free domains 11(77) and n(J r ) for the parameters 
s and t, with u = /3. Right: Detailed view of n(J r ), showing the trapezoidal domain 
II*(J r ) of equation (211. The latter is constructed as the intersection of II (Jr) (semi¬ 
transparent green trapezoid), n(7 r 2) (semi-transparent blue triangle), and n(J r 3 ) 
(semi-transparent red square). The seven periodic tiles Vi do not contribute any 
additional constraints. The unperturbed one-parameter model corresponds to the 
(yellow) line s — t = /3 2 , which lies along the south-east boundary of n(J r ) 


boundaries of merged atoms), plus a finite number of periodic tiles. 

For the plane u = 0 of the five-atom map, the case N = 1 appears to be 
the norm, so that in those models the exceptional set is likely to be uniquely 
ergodic. Elsewhere, a proliferation of ergodic components is typical (albeit 
not universal). Notably, we have found no example of a rigidly self-similar, 
single component piecewise isometry with two or more parameters. Whether 
such a dynamical system exists at all remains an important open question. 
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